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Experiments on dry granular matter flowing down an inclined plane are performed in order to 
study the dynamics of dense pyroclastic flows. The plane is rough, and always wider than the flow, 
focusing this study on the case of laterally unconfined (free boundary) flows. We found that several 
flow regimes exist depending on the input flux and on the inclination of the plane. Each flow regime 
corresponds to a particular morphology of the associated deposit. In one of these regimes, the flow 
reaches a steady state, and the deposit exhibits a levee/channel morphology similar to those observed 
on small pyroclastic flow deposits. The levees result from the combination between lateral static 
zones on each border of the flow and the drainage of the central part of the flow after the supply 
stops. Particle segregation features are created during the flow, corresponding to those observed on 
the deposits of pyroclastic flows. Moreover, the measurements of the deposit morphology (thickness 
of the channel, height of the levees, width of the deposit) are quantitatively related to parameters 
of the dynamics of the flow (flux, velocity, height of the flow), leading to a way of studying the flow 
dynamics from only measurements of the deposit. Some attempts to make extensions to natural 
cases are discussed through experiments introducing the polydispersity of the particle sizes and the 
particle segregation process. 
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INTRODUCTION 



Flows of grains are encountered in many natural sit- 
uations: sand dunes, debris flows, landslides, turbidity 
currents, avalanches... These flows involve complex in- 
teractions between grains, and sometimes a fluid. For 
example, the dynamics of pyroclastic flows can follow 
every possible combination between two end-members, 
depending on relative intensities of particle-particle, and 
gas-particle interactions: dilute suspensions (only gas- 
particle interactions), and granular flows (only particle- 
particle interactions). If large volume pyroclastic flows 
are related to gravity currents (dense or dilute) 0, |E 
3], small volume flows behave as granular flows, pro- 
duced either by dome or by eruptive column collapse 
IE IE IE IE HI, P Many models have been proposed 
(reviews in [HJ, El ) , ranging from dilute turbulent sus- 
pension flows dominated by fluid-particle interactions 
El IE 113 EE EE EE- to dense gravity currents (con- 
centrated suspensions) with fluidization effects and some 
particle-particle interactions [iH ITtL ITsl ll^. Eol l2l | . and 
to flows controlled by particle collisions [21|. But, to our 
knowledge, there is no study involving frictional interac- 
tions between particles. 
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In this paper, we focus on dense dry granular flows. 
They are defined as a flowing layer of grains, in which 
air-grain interactions are negligible (consequently, heat 
exchanges do not influence processes). For dense granu- 
lar flows, interactions between grains are dominated by 
friction |22ll23| and not by collisions, as it is the case for 
rapid granular flows |24l Y2^L 1261] . Moreover, if the flow is 
said "dry" all cohesive forces between the grains are as- 
sumed negligible. Their dynamics is different from those 
of gravity currents which obey fluid mechanics, and in 
which grains are, at least partly, suspended by the ver- 
tical component of the fluid velocity. Gravity currents 
loose particles by sedimentation, and get dilute, unlike 
dense granular flows whose volume fraction of grains re- 
mains close to 60% . Some pyroclastic flows have 
intermediate dynamics (a granular layer flow below a 
gravity current), and pyroclastic granular flows are sur- 
rounded by a suspension of ash. The aim of this paper 
is to obtain laws for the dynamics of the 'dense granular 
flows' end-member, and quantitative relations between 
their dynamics and the morphological characteristics of 
their deposits. In that way, it should be possible to get 
quantitative information about flow velocities from pyro- 
clastic deposit measurements. 

The assumption about the occurence of dense pyro- 
clastic granular flows, or their dense granular basal part, 
is mainly supported by deposit observations: segregation 
features (fig. QJ, erosion due to the friction between clasts 
|27| , and erosion of the subtratum due to the dense char- 
acter of the flow 0,H^|. Moreover, there are no pyroclas- 
tic deposits on high slopes although they appear when 
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slope decreases 0, |(| |2^. When not laterally confined 
in valleys, the deposits display a particular morphology: 
the lobes have high parallel lateral levees (about 1 m 
high) enriched in large blocks while the central chan- 
nel is lower and mainly composed of smaller particles 

IllllllllllillillflEllalESISlEl The fronts 

have rounded bulbous shapes and steep margins (fig. ^) 
0, l28l l33| . Lobes deposited on the lowest slopes exhibit 
similar shapes. Various hypothesis have been proposed 
to explain the levee/channel morphologv: lateral static 
zones linked to a Bingham rheology [|| |2^, |^| , drainage 
of the central part of the deposit with static levees [33| , 
differential deflation after emplacement |33j and differ- 
ential fluidization of borders during emplacement 
both due to the heterogeneity in the spatial distribution 
of particle sizes. This includes interaction of the borders 
with ambiant air creating differential flotation and seg- 
regation |23, and rafts of pumice clasts stranded on the 
sides of the flow 0. One aim of this paper is to know 
which part of these morphology and segregation features 
is due to the dense flow of dry granular matter. 




FIG. 1: Lobe-shaped deposits produced by the 1993 pyro- 
clastic flows of Lascar volcano (Chile): (top) levees enriched 
in large pumices (flow moving ahead), and (bottom) global 
view of lobes observed on the South-East flank of the vol- 
cano. The cube gives a 3D scale: edges are 78 cm long (photo 
from authors). 

Most experimental studies of dry granular flows con- 
cern confined flows: these flows are either confined in 
narrow channels [2(| |3(| 153, IH, |3j|, or in a wide ge- 



ometry set-up where friction at lateral walls is negligible 
|4fil | . The width of the flow is imposed by the width of the 
experimental set-up: the flow and deposit morphologies 
are strongly affected by confinement (in a similar way, ex- 
cept a few studies on snow avalanches |41| . experiments 
on gravity currents have been done in channels, in which 
the transverse morphology of the deposit is impossible to 
study). Some studies focus on the spreadin g of a finite 
mass of granular matter on a wide slope |42ll43j] , but no 
steady state regime is reached, and the deposit morphol- 
ogy is strongly influenced by its initial shape, or deformed 
by local changes of the slope [44ij . Experiments involving 
the destabilisation of a bump of material, irregularly fed 
during about a quarter of the flow duration, show that 
the levee/channel morphology can develop [43. But, to 
our knowledge there is no data on dynamics of uncon- 
fined dry granular flows on inclines, that are fed with 
a constant flux for a long duration. Models based on 
kinetic theory HE E2| are valid for rapid collisional 
flows, although depth-averaging equations with a basal 
friction seem more appropriate for dense frictional flows 
[4fT | (review in [2]|). All these studies show that there 
is no layer with a zero vertical gradient of the horizontal 
velocity at the surface of the flow, eliminating a Bingham 
rheology. 

In this study, laboratory experiments involving dense 
granular flows have been performed in order to relate pos- 
sible flow regimes and associated deposit morphologies. 
The input flux is constant such that the flow reaches a 
steady state. The flow is never as wide as the set-up, and 
its width is free to adapt to the dynamical conditions. 
Precise experiments have been carried out to understand 
the mechanism of formation of the levee/channel mor- 
phology and to quantitatively link information on the 
flow dynamics and measurements of the deposit morphol- 
ogy. 



II. EXPERIMENTAL APPARATUS 

The flows are created by release of microbeads on an in- 
clined plane (fig. [5J . The plane (80 cm wide, 2 m long) is 
made rough by gluing 425-600 /im beads on it. The slope 
is chosen such that the flow is dominated by frictional 
interactions between the beads, and between the beads 
and the plane, and not by collisions [25ll46| . Two feeding 
systems are used: (1) A funnel with a variable aperture 
allows control of the flux imposed by a point source pro- 
ducing unconfined finger-shaped flows. (2) A box with a 
40 cm gate is used to create sheet flows, rapidly becoming 
laterally confined, i.e. 80 cm wide flows. 

Particles are 2.5 g/cm 3 density glass beads. Their di- 
ameter is either in a narrow (but not monodisperse), or in 
a wider (named 'polydisperse' here) range of sizes. The 
humidity in the laboratory is held between 50% and 55% 
to avoid electrostatic effects or capillary bridges between 
particles, which are cohesive effects |47| . 

The morphologies of flows and associated deposits are 
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FIG. 2: (top, left) A glass beads finger-shaped flow on a rough 
80 cm wide and 2 m long inclined plane, and (top, right) 
the associated deposit lighted from the right, showing lateral 
levees and low central channel. The width is around 12 cm 
in each case, (bottom) The deviation of the laser sheets, pro- 
jected perpendicularly to the flow direction, is proportional 
to the deposit and flow thicknesses. The scale is given by the 
calibrated holes on the small plate (right of the picture); the 
width of this deposit is around 30 cm. 



measured by projecting several laser sheets on the plane, 
perpendicularly to the flow direction (fig. |3 bottom). 
The deviation of each line is proportional to the thickness 
of granular material on the plane. A reference allows the 
scaling of the transverse profiles. 

For each set of experimental conditions, the width L 
and the thickness h of the flow as well as the deposit 
morphology are measured. The propagation of the flow 
is followed in time, giving the front velocity v, that cor- 
responds to the mean flow velocity. 



III. THE LEVEE/CHANNEL MORPHOLOGY 
AND OTHERS 



First, the mode of formation of levees in granular flow 
deposits is presented. Then, the observed flow regimes 
are described and related to the morphology of the as- 
sociated deposits. Finally, we show how it is possible to 
get quantitative information on the flow dynamics, using 
only measurements of the deposit morphology. 



A. Formation of levees and channel 

Particles constituting the flows are in a narrow range of 
sizes (typically 300-400 //m) . The flow develops from the 
point source. When moving downwards, it spreads lat- 
erally until it reaches a constant width. Once its stable 
width reached, the position of the front varies linearly 
with time: the velocity is constant (fig. |3J). The flux 
being constant, the thickness of flowing material is also 
constant and the flow reaches a steady state regime (par- 
tial derivatives of velocity, width, height, with respect 
to time are equal to 0). When fully developed, the flow 
is finger-shaped with parallel borders and a curved bul- 
bous front (fig. |2 top, left), corresponding to the shape 
of deposits in the field [34|. 
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FIG. 3: For fluxes of: (A, A) 6 g/s, (♦, 0) 8.5 g/s, (•, O) 19 
g/s, and (■. □) 34 g/s: (plain symbols) longitudinal variation 
of the width of the flow down the slope (length is the distance 
from the source) . The flow reaches its stable configuration at 
a distance between 50 and 150 cm from the source depending 
on the flux; (open symbols) variation of the position of the 
front of the flow (i.e. length of the flow) versus time. The 
linear evolution, after a first slightly curved part, shows that 
the velocity of the flow is constant once it has reached its 
stable width (beads: 350-390 /im, roughness: 425-600 /im, 
slope: 25°). 

The evolution with time of the flow morphology has 
been followed at a distance chosen such that the flow has 
reached its stable configuration. A 'small' flux (see sec- 
tion is used for technical reasons. Three distinct 
phases occur (fig. 0J: 

- First, the front of the flow arrives at the chosen distance: 
the height and the width of the cross-section increase un- 
til a stable profile is reached. 

- Second, the profile does not evolve. After a small de- 
crease of the height corresponding to a slight longitudinal 
convex shape of the front, width and height are constant 
with time, at the chosen distance. They are also con- 
stant along the flow direction. We observe lateral static 
zones on each border of the flow: pictures of the laser 
lines taken with a long exposure time are blurred where 
beads are moving, and not blurred where there is no sur- 
face motion (limits drawn in fig. |3J. The static zones 
can also be seen when a line of dark particles is poured 
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instantaneously on the surface of the flow (inset). The 
line, initially perpendicular to the flow direction, evolves 
showing the surface velocity profile: the dark line is de- 
viated neither on the plane, nor on the lateral borders of 
the flow. 



A small volume of colored particles was also suddenly 
released from the funnel and their trajectories were 
followed on the flow surface. Roughly, particles reaching 
the center of the front are re-injected into its rear, 
while particles reaching the lateral borders of the front 
become static, building the levees (fig. [5J|. Thus, the 
composition of the static zones records, down along the 
slope, the temporal evolution of the composition of the 
front. A thin line of colored particles was also placed on 
the plane, on the future path of the flow. These particles 
reappeared on the surface of the flow, at the rear of the 
front, showing that the recycling process involves the 
whole thickness of the flow, and that the basal part of 
the flow is not static. 

- Finally, as the supply stops, the height in the channel 
decreases due to the downward flow (fig. while 
the height of the static borders remains unchanged. 
Consequently, the central channel becomes lower than 
the borders, leading to the levee/channel morphology. 
More accurately, the static borders partly collapse into 
the channel when it empties. The summits of the 
levees are slightly less high than the maximum height 
of static zones, and also a little further apart than the 
limits of static zones. The levees record the maximum 
height reached by the flow at the future location of their 
summits. They are thus linked to the height of lateral 
static zones just behind the front, and consequently 
linked to the front height. Here, we can associate the 
height of levees with the height of the flow, because 
the longitudinal convex shape of the front is small. 
But, when the front has a stronger convexity (case of 
some polydisperse media), it would be more accurate to 
associate the height of levees and the maximal thickness 
of the front. 



Consequently, the formation of levees is associated 
with the existence of lateral static zones and the partial 
emptying of the central portion of the flow. This is com- 
patible with the precise observations and interpretation 
given by [33J. 

When supply stops, the flow continues to propagate on 
a short distance, using the flux coming from the partial 
emptying of the back of the deposit. The deposit ex- 
hibits a smooth transition between a rounded front and 
the levee/channel morphology at the back, of which the 
deposit mainly consists (fig. right). This morphol- 
ogy is comparable to those observed in the field |23. It 
has also been experimentally observed in the failure of a 
bump of sand |45|. 
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FIG. 4: (top) Arrival of the front of the flow at the chosen 
location: both height and width increase with time. X-axis is 
taken perpendicularly to the flow direction, (middle) Stable 
flow with constant width. The slight decrease of the height 
corresponds to a longitudinal convex shape of the front. The 
two lines limit the area where there is a surface motion: lat- 
eral static zones exist on each border of the flow (slightly 
wider at 18 s). Inset: a line of dark particles poured instanta- 
neously on the surface of the flow shows the two static lateral 
borders and the surface velocity profile, (bottom) After stop- 
ping the supply, the central channel partly empties: the final 
morphology is high lateral levees with a low central channel. 
Here, the channel is not flat because of the 'small' flux used. 
The levees record the height reached by the flow in the lat- 
eral static zones (beads: 300-400 /im, roughness: 425-600 (jm, 
flux: 7 g/s, slope: 25°) 
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FIG. 5: Schematic representation of particle trajectories at 
the flow surface: particles arriving at the front, either go to 
build the lateral static borders at the rear of the front, or 
are recycled in the front when close to the center of the flow. 
Consequently, the levees record, down along the slope, the 
temporal evolution of the front composition. 




FIG. 6: (left) A local variation of the slope from 23.2° to 
29°, on a mean slope of 26°, induces a local enhancement of 
the morphology; (right) after the supply stops, the front of 
deposit exhibits a rounded shape turning progressively into 
the levee/channel morphology at the back. 



B. Evolution of the levee/channel morphology 

1. Evolution with the angle of the slope 

The morphological aspect ratio varies when the incli- 
nation of the plane 9 changes (fig. 00 top). Because both 
levees height, hi evi s e , and channel thickness, h c \ lanne i de- 
crease with 9, h, 
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I being inferior to /i/euee, the mor- 
phology is enhanced for steep slopes, as observed in the 
field |3l| . Simultaneously, width and velocity increase 
with 9 (fig. 0). 

The dynamics of confined dense granular flows down 
rough inclines show that in a range of inclinations (6±, 
$2), a uniform layer of deposit remains after a flow has 
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FIG. 7: (o) width of flow and deposit, and (♦) mean flow 
velocity (i.e. velocity of the front) both increase with the angle 
of the slope (flux: 7.1 g/s, beads: 300-400 fim, roughness: 
425-600 fim). 



passed ^3,H^|- The deposit thickness, named h stop , is 
independent of the velocity and of the thickness of the 
flow h, but depends on 9 (fig. EJ). The thickness h stop 
is interpreted as a characteristic length on which stresses 
generated by the friction with the substratum are trans- 
mitted in the granular matter. It is linked to the precise 
characteristics of the plane roughness and of the particle 
size . h s top is one possible vertical length scale for the 
system |40| . 

Previous authors found that there is also a maximum 
thickness of material which can stay on an inclined plane 
without moving, h start 0, 0] ; h s t ar t depends on 9, 
generating a similar curve, shifted towards large angles 
{hstart > h s top)- A layer of granular matter with a thick- 
ness smaller than h s top is stable and never flows. If its 
thickness (or 9) increases and reaches h s tart(9), the layer 
destabilises and flows; the produced deposit having the 



hstop thickness. 



For a thickness between h s top and h s 



the granular layer can be stable or unstable depending 
on the initial conditions (hysteresis). 

The solid line in the figure El is a fit for h s top{9) estab- 
lished by Pouliquen i 4($ where C is an adaptable param- 
eter and d the size of the beads: 



tan 9 = tan 0i + (tan 9-2 — tan#i) exp 



hstop 
Cd , 



(1) 



In our experiments, h c h a nnei and hi ev ^ e show a de- 
creasing trend compatible with a l hstop{9)' type curve, 
although hi eV e e could be an i h s to P {9) + constant' type 
curve (fig. EJ- However, we do not have enough data 
to check that point. This shows that the thickness of 
deposit in the channel is linked to the stopping process, 
and is independent of height and velocity of the flow. 

To confirm that h c i ia nnei corresponds to the h s top pa- 
rameter defined by Pouliquen 40], we measured the 
thickness of our material lying on our plane after a wide 
confined flow has passed. The flows were made using the 
40 cm gate, sufficiently opened to obtain a thick sheet 
flow, which rapidly reaches the edges of the plane, giving 
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FIG. 8: (top) Increasing the angle of the slope enhances the 
morphology of the deposit. From top to bottom: 25°, 27°, and 
29° (flux: 7.1 g/s). (bottom) Evolution of the morphology of 
a flow and its deposit for different fluxes, from top to bottom: 
3.7, 8, 17.5, and 38 g/s (slope: 25°). The width of the levees 
varies slightly, although the width of the deposit varies a lot 
(beads: 300-400 fim, roughness: 425-600 fim). 



a confined flow and a flat deposit. The measurements 
are less precise because the scale cannot be on the same 
picture as the deposit. Nevertheless, data on both thick- 
nesses are in good agreement (fig. EJ- 

In the field, the levee/channel morphology has been 
described by the ratio h c hannei/hi eve - e . This ratio is sen- 
sitive to errors: we fit both sets of data, and calculate the 
ratio of each fit (fig. 0). The decrease of the ratio with 8 
is compatible with observations on pyroclastic deposits: 
a ratio of 0.9 for 'low' slopes, and of 0.5 for 'steep' slopes 
on 1980 Mount St. Helens deposits |3J]. Moreover, the 
morphology adapts very rapidly to a local variation of 
the slope, leading to a local deformation of the deposit 
(fig. El left), which has been observed in the field 
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FIG. 9: Thickness of deposit: (O) hi ev &e, and (■) h c h a nn<d- 
Fits come from equation Q with the same values of 8i and 82 
according to: ( — ) "h s to P curve", and (....) u h s t op + constant" 
curve. (♦) the h c h e nnai / hievee ratio and (- -) the fit obtained 
from the ratio of both fits show the decrease observed in the 
field. ( X ) Thickness of the deposit of wide laterally confined 
flows, with error bars. Thicknesses of the deposit of both an 
unconfined flow, and a confined flow are equal (flux: 7 g/s for 
confined flows, and large enough to cover the whole plane for 
confined flows, beads: 300-400 fim, roughness: 425-600 |im). 



cept for very narrow flows, the morphology of flow and 
deposit is mainly flat. This is appropriate to define h 
as the thickness of flow and h c hannei as the thickness 
of deposit in the channel. For narrow flows (fig. |3J , h is 
taken as the maximum thickness of the flow, and hchannei 
as the minimum deposit thickness. Second, the width of 
the levees is almost the same, independently of the width 
of the flow, i.e. independently of the flux. Consequently, 
the aspect ratio changes as the flux varies. This morpho- 
logical change is different from those induced by a slope 
variation (fig. [S] top) , giving a way of differentiating the 
effects of flux and of topography in the field. 



C. Flow regimes and associated deposit 
morphologies 

The second step of this study consists in linking the 
flow dynamics and the deposit morphology. A qualita- 
tive link is presented first, then a quantitative relation is 
determined for one regime and, finally, the influence of 
polydispersity of materials is discussed. 



2. Evolution of the morphology with the flux 

A flux increase induces an increase of the flow veloc- 
ity and width (fig. [3}. At the same time, the thickness 
adapts. Figure |S] (bottom) shows a comparison of flow 
and deposit cross-sections for different fluxes. First, ex- 



1. Description 0} the regimes 

This study is focused on dense frictional flows. The 
collisional regime has not been studied, although it is 
observed for higher slopes. Several flow regimes happen 
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depending on flux and inclination of the plane. The 
volume fraction of beads is between a random loose 
and a random close packing of beads (0.56-0.64%) |Hcj . 
and if there is any difference between the regimes, they 
would be to small to be measured. Each type of flow 
has its own dynamics, and produces a particular type of 
deposit (figs. 121 and ITU)l . Three regimes can be defined, 
depending on values of slope and flux (fig. In this 

section, experiments involve particles with narrow size 
distributions (215-224, 150-250 or 300-400 fim), on 425- 
600 /im roughness. This has been done to obtain various 
values of the flux. Results are consistent even if beads of 
different sizes are used on the same substratum, because 
the size changes only slightly. A few test experiments 
have been done with rough painted beads, and with 
other bead sizes. They show a dependency of the limits 
of the regime fields on the characteristics of the particles. 





FIG. 10: Two other types of (top) flow and (bottom) asso- 
ciated deposit: (left) intermittent flow regime and rounded 
shape deposit (for low flux and/or slope), and (right) roll- 
wave regime and thin flat deposit (for high flux and slope). 
The line is the projection of the laser sheet either perpendic- 
ular or parallel to the flow direction. 



and section UlI Afl . The flow is finger-shaped, with con- 
stant width, height, and velocity. The deposit displays 
the levee/channel morphology. For this regime, a 'small' 
flux corresponds to conditions close to the limit between 
intermittent and steady state regimes. 
-Finally, for high flux and slope, a roll-wave regime oc- 
curs (figs. 1101 right and II 1|) . The flows are thin and 
finger-shaped with parallel lateral borders which are not 
very-well defined. Small waves develop at the surface of 
the flow and rapidly reach its front. The deposit is cither 
flat and as thin as 1-2 layers of beads, or non-existent. 
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FIG. 11: Flow regime fields: (□) intermittent flow regime (see 
fig. 1101 left), (+) roll- wave flow regime (see fig. 1101 right), 
and (•) steady state flow regime (see fig. |2J. The curve comes 
from the h 3 t op (8) curve (( — ) in fig. [5J|, adapted with a linear 
relation between h s to V and the flux. 



-First, for low flux and/or low slope, an intermittent 
flow regime is observed (figs. left and lllfl . A rounded 
shaped flow spreads down the slope, then stops. As it is 
continuously supplied with granular matter at the back, 
a second flow starts a few seconds later from the accu- 
mulation zone, and covers the previous one, propagating 
the whole deposit a little further, and then it stops again. 
In the figure El (top, left), the front is static, while the 
back is moving ahead (blurred area). The flow is con- 
sequently intermittent, with a front velocity alternating 
between and a positive value, which might be similar to 
features observed during Montserrat eruptions |30|. No 
measurement of frequency has been done, which seems 
to depend on both flux and slope. When the deposit is 
long enough, several successive avalanches can propagate 
simultaneously. The deposit has a rounded shape (fig. 
1101 bottom, left). Sometimes, when the feeding has been 
stopped early, it exhibits a step, the provided volume be- 
ing not large enough for a total covering. 
-Second, for a given intermediate range of flux and in- 
clination, a steady state flow is observed (figs. [21 [HI 



These three regimes are closely related to those ob- 
served for wide confined flows |40l l5l| , when making an 
equivalence between the flux and the flow thickness: for 
low angle and/or thickness, no flow occurs; for interme- 
diate values, flows are in a steady state regime; for high 
angle and thickness, instabilities appear, which are pre- 
dicted by the theory Using sand, these waves 
appear as soon as there is a flow [HlJ , and could be asso- 
ciated with the levee/channel morphology phf . 

These regimes obey different dynamics. It is important 
to know in which regime the flow responsible for a deposit 
has propagated, before any attempt to calculate its ve- 
locity. Fortunately, the deposits have different morpholo- 
gies, and the flow regime can be deduced from qualitative 
observations of the deposit. Note that the determination 
of the regime is clear when the levee/channel morphol- 
ogy is observed. But when the deposit is rounded, it 
can be either the result of an intermittent regime, or the 
result of the propagation of the flow within the steady 
state regime, the observation being made too close to 
the front to see the emptying of the back of the deposit 
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(fig. El right). 



2. Interpretation using the thickness h a top{9) 

It is possible to interpret the regimes using the thick- 
ness hstopiO) remaining after a wide confined flow has 
passed (fig. and eq. [TJ. Two angles are involved: 
for 9±, h st0 p tends towards infinity, while above 02 , it 
is equal to zero. This means that below 9i, no uniform 
layer of granular matter can be thick enough to flow, and 
that above 82, no deposit remains (even if some matter 
may rest on a slope slightly above 62, when placed in a 
static state on the plane). The dashed curve (fig. 
comes from the h s t op (0) of the figure El assuming that 
the flux and h stop are linked by a linear relation. That 
shows that the limit between steady state and intermit- 
tent regimes, and the h s t op (9) curve both go through a 
compatible range of angles [9\, 62) ■ 

-Under the h s ± ov (9) curve, a wide confined granular layer 
does not flow [4£j . When granular matter is released from 
a funnel, it cannot initiate a flowing layer with a local 
constant thickness, either because the slope is below 9\, 
or because the flux is too low to feed a flow thicker than 
h s top- K constitutes a heap till its free surface reaches an 
angle above the stability angle. Then, the heap collapses 
and the granular matter flows in a non-symmetrical way 
on the slope. The intermittent regime is the result of 
successive avalanches building a 'pile' on a steep slope, 
because a uniform layer would not flow. 
-Above the h stop (9) curve, the amount of released matter 
is sufficient to feed a layer which is locally thicker than 
h stop : thus the steady state regime occurs. 
-For angles slightly above 6*2, the deposit consists of two 
levees aside an empty channel, because for every angle, 
the thickness of the levees is larger than those of the 
channel (fig. EJ- 

-Above 6*2, the particles flow down and no deposit stays 
on the plane. Extrapolating our results, it seems that 
both the steady state and roll-wave regime can occur. It 
is also very likely that the collisional regime occurs. 

Consequently, the potential occurrence of the different 
types of deposits can be located in comparison with the 
h s top{0) curve (fig. [T2J) . 

-For slopes smaller than #1, the deposit is thick and 
rounded in shape. 

-For slopes between 9\ and 6*2, the three regimes can oc- 
cur. The deposit displays the levee/channel morphology 
except close to the front, or except if the flux is too low 
or too high. Very close to #2, only a few particles stay 
on the plane (h s t op < d), corresponding to lag deposits 
observed on steep slopes [3^ |. 

-For slopes greater than 9 2 , there is no deposit (at least 
in the channel), as it is observed in the field 0,1a, u3- For 
angles slightly above 6*2, the "levees with empty channel" 
morphology exists (fig. El right inset). 

We conclude from these results that if the 
levee/channel morphology is observed on a deposit, 
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angle of the slope 



FIG. 12: Location of the deposit types, shown by the inset 
photographs, compared to the h s to P (6) curve. For flow condi- 
tions below the curve, the deposit has a rounded shape. The 
levee/channel morphology exists above the curve except when 
instabilities occur (high angle and thickness). 'Levees with an 
empty channel' can be observed for angles slightly above 02 
(dark beads for this inset). For higher angles, there is no de- 
posit either in the levees, or in the channel. Consequently, 
the levee/channel morphology, with some matter remaining 
in the channel, occurs only between the 9\ and 82 angles. 

there is no doubt that the associated flow was in a 
steady state frictional regime, and that the angle of the 
slope was between 9\ and 9 2 - We made observations on 
the south-east lobes of the Lascar (Chile) 1993 deposits, 
which gives approximate values for 9\ and 02 of 6° and 
14°. These values correspond to erosion/deposition 
behaviors observed for slopes respectively equal to 
14-16° and 4-6° at Montserrat 30]. Both angles are 
small, probably due to the large values of the ratio 
between substratum roughness and particle sizes |49j . 

Note that, for a given erupted volume, knowing 
h s top{9) is not enough to calculate the longitudinal de- 
posit extension (except, if the flow in confined in a val- 
ley): it could be long and narrow (small flux) or short 
and wide (large flux). 



D. From deposit morphology to quantitative 
information on flow dynamics 

There is a good analogy between confined and un- 
confined flows, confirmed by the existence of analogous 
regimes [HlJ and we are confident to try successful meth- 
ods used for dynamics of confined flows for the interpreta- 
tion of dynamics of unconfined flows. The challenge is to 
obtain some quantitative values for the dynamics (flux, 
velocity, height of the flow) only from measurements per- 
formed on the deposit. 

We focus on the levee/channel morphology because the 
particular shape of the deposit can be related without 
ambiguity to the flow regime. In this regime, the influ- 
ence of the flux has been tested. It appears that the 
larger the flux, the faster the front of the flow and the 
wider the flow (and the deposit) (fig. I13fl . As expected 
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|40t l4Sj . hchannei contains no information on dynamics, 
which is not the case for hi eve - e and h (fig. I13fl . Note that 
the variation of hi eV ee with the flux is incompatible with a 
2-dimcnsional calculation of the dynamics of an infinite 
layer of fluid having a Bingham rheology, which would 
predict that the thickness of the static zone is determined 
by the yield stress and should depend only on the slope. 
The levee height is linked to h, through the formation of 
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FIG. 13: (top) (o) Width of flow and deposit, and (♦) mean 
flow velocity (i.e. velocity of the front) both increase with the 
flux, (bottom) Thickness of (•) flow and (■) levees increases 
slightly with the flux in a similar way. Error bars come from 
the difference between both levees heights in each deposit, (o) 
The thickness of material in the channel of the deposit almost 
does not depend on the flux (beads: 300-400 ^m, roughness: 
425-600 fim, slope: 25°). 

lateral static zones, and to h c hanneU through the partial 
collapse. Because the h c hannei thickness is constant (or 
very slightly decreasing) with the flux, the increase of h 
and hi eV ee with the flux are similar. More accurately, the 
slightly larger value of hchannei for small flux may be due 
to an arching process between the static borders during 
the collapse |53t l55j . Because the associated flow is 
narrow, the short distance between the static zones pre- 
vents the complete emptying until the h s t op thickness of 
a wide flow deposit is reached. It would have been more 
accurate to keep the notation l h s t p' for deposit produced 
by a wide flow, and to say that the thickness in the chan- 
nel, h c hannei tends to h stop when L tends towards infinity. 



we chose not to distinguish between these 2 variables. 
The morphology of the deposit then contains parameters 
linked to the flux (hi eV ee , L) and parameters independent 
of the dynamics (h s t op ). 

For wide laterally confined flows down rough inclines, 
Pouliquen ^(j nrs t established a relation between the 
mean flow velocity v and the flow thickness h (eq. [21 ■ As 
previously, h s t p is the thickness of the uniform deposited 
layer. It is the vertical length scale of the system and used 
to normalize the flow thickness: 



gh 



^stop 



(2) 



g is the gravity, and (3 a constant equal to 0.136 for glass 
beads flowing down a plane made rough by gluing other 
glass beads. Due to the normalization by h s t op , all flows 
composed of different sizes of glass beads, on different 
planes of this type, follow the same law. This type of 
correlation is also verified in our experiments on uncon- 
fined flows, giving a similar value for /?. 

Figure El shows the correlation obtained in the case 
of our experimental flows by using hi eV ee instead of h, 
and h c hannei for hgtop- As expected, the value of the 
constant f3 is different from those previously obtained, 
being equal to /3—Q.5 (fig. I14JI . The correlation exists, 
although the slope is high, showing that the correlation is 
of the same order of magnitude than the data variability. 
Unfortunately, due to the finite length of the experimen- 
tal set-up, it was not possible to study the steady state 
regime for higher fluxes, and consequently to accurately 
determine this correlation. Nevertheless, we deduce a re- 
lation allowing an estimate of the flow velocity only using 
deposit morphology: 



Vghie 
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FIG. 14: Normalized correlation between velocity and thick- 
nesses of the levees and the channel, for laterally unconfined 
granular flows and associated deposits (beads: 300-400 /im, 
roughness: 425-600 /im, slope: 25°). 
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Its extension to natural systems is not direct since the 
value of the constant (3 for volcanic granular material is 
unknown and probably different (experiments using sand 
give /3=0.65 |51|). We are also aware of the deflation 
of the deposits, which is low in most cases (0.3-1% for 
St Helens deposits j^H or 2-6% for Unzen deposits 
but can also be important (30-50% for St Helens deposits 
|33j ) and from which the deposit measurements should be 
corrected before using this relation. Each deposit have 
to be corrected in a different way and it is not possible to 
obtain a general equation including the deflation effect. 
The deflation is not calculated in this paper, since, to 
our knowledge, no precise morphological natural data are 
available. Moreover, natural materials have a wide size 
distribution inducing segregation processes that are not 
taken into account in relation^ and whose influence will 
be shown in the next section. 



E. Flows with a wide particle size distribution 

During experiments involving wider size distributions, 
segregation features have been observed on flow and de- 
posit. The large particles rapidly reach the surface of the 
flow and consequently accumulate at the front and in the 
static borders. In the deposit, the levees and the front 
are enriched in large particles while the channel contains 
mainly small particles, as observed in the field (figs. ^ 
d]and 0)- 

Segregation is also responsible for fingering of the flows 
observed in the laboratory as well as in the field (black 
arrow in fig. 1151 bottom). The jamming of a large par- 
ticle when it reaches the substratum at the front of the 
flow induces the flow to separate into two distinct lobes 

na. 

The influence of the size distribution on the velocity 
of the front and on the morphology of the deposit was 
also tested. It results from two processes: (1) contacts 
between particles in the flow are modified, i.e. the rhe- 
ology changes; (2) the basal friction is strongly modified 
depending on which particle size is in contact with the 
plane. The friction depends on the size ratio between 
particles and substratum roughness. For monodisperse 
roughness and grain sizes, it has been shown that the 
friction is smaller for high and small ratios, than for ra- 
tios close to 1 49]. For that reason, h in equation (J5J) 
needs a normalization by h stop . Because of segregation, 
fine particles accumulate in the basal part of the flow, 
and eventually, fill the holes, smoothing the roughness. 
The plane thus becomes smoother for a large particle, 
but could be more rough for a small one. The friction 
between flow and substratum is mainly the result of the 
interaction between these small beads and the plane and 
could be reduced or enhanced depending on which sys- 
tem it is compared to (pure flow of small particles or 
pure flow of large ones), and depending on the ratios of 
particles and roughness sizes. 




FIG. 15: (top, left) Accumulation of large black particles in 
the front of laboratory granular flows, and (top, right) accu- 
mulation of large pumices in the front of the right lobe, in 
the Lascar pyroclastic flows deposits. Lobes fingering (bot- 
tom, left) in a glass beads flow deposit and (bottom, right) 
in a Lascar pyroclastic flow deposit (main flow going to the 
far right). The black arrow shows the point where the flow 
separates in two lobes. In the left picture, the small flow is 
particularly rounded in shape and still covered by large par- 
ticles, because it did not propagate far enough to empty the 
back of the deposit. There is an accumulation of large black 
particles and of large pumices in the levees (see also fig. Q 
top). 



1. Bidisperse media 

The influence of fine particles has been observed in 
bidisperse flows involving a small fraction (respectively 5 
and 10 weight%) of 45-90 fim beads in 300-400 fim beads. 
For such sizes, the fines fill the holes of the plane consti- 
tuted of 425-600 ^m glued beads. The velocity increases 
compared to the flow of pure large particles (by a fac- 
tor of, respectively, 1.4 and 1.9), and the morphology is 
strongly enhanced: the channel is almost empty between 
two high levees (see fig. right inset). 

On the same plane, each particle size is associated 
with a different quasi-parallel h sto J9) curve |4^. F° r 
these sizes and roughness, a model ]4!j predicts that the 
h s top(9) curve of the fines is shifted to small angles com- 
pared to the h st op(0) curve of the large ones. In bidisperse 
mixtures, segregation puts the fine particles at the bot- 
tom of the flow, and the thickness of the deposit in the 
channel is mainly imposed by the h s t op {0) curve of the 
fine particles. But, pure large particle flow and mixture 
particle flow would roughly lead to similar levees height, 
levees being built at the rear of the front, mainly con- 
stituted of large particles in both cases. At the present 
angle (27°), small particles fill in the channel during the 
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flow, but the channel deposit is almost empty, because 
the slope is close to 82 of the fines particles. On the 
contrary, the channel contains some matter for a flow of 
pure large particles, the present angle being far under 82 
of the large particles. Moreover, a small value of h stop 
corresponds to a low basal friction: the velocity of the 
bidisperse flow is higher than the velocity of a pure large 
particle flow. The shift of the h s t op (8) curve for the small 
particles is responsible for enlarging the range of angles 
where the "levees with empty channel" morphology oc- 
curs. This morphology has been observed at least from 
25° to 30°, and the associated levees are not necessarily 
small, which is not the case of monodisperse flows. 
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2. Continuous polydisperse media 

Continuous polydisperse granular media have also been 
used. Successive mixtures (table I) have been made, cor- 
responding to a progressive increase in the width of the 
size distribution. 



TABLE I: Composition of the polydisperse mixtures: abun- 
dances are in weight% (left column: size ranges in fim; top 
line: batch number used in the figures). 
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Morphologies of flows and deposits are enhanced when 
the size distribution is wider (fig. [To]l . This trend quickly 
appears and is linked to the segregation process combined 
with the shift of h s top curves when changing the smallest 
and the largest sizes in the distribution. Observations 
on deposits show that when the size distribution is in- 
creased, the particles of the channel become smaller and 
the particles of the levees larger. 

Consequently, a variation of the height of the levees, 
or of the thickness of the channel can be explained either 
by a variation of the flow composition (fig. I17|) or by a 
variation of the flow dynamics (fig. I13|) . The normalized 
relation between velocity and deposit morphology (eq. 
has to be extended for polydisperse media before being 
used to quantitatively relate the dynamics and morphol- 
ogy in the field. 

On the other hand, the velocity increases only very 
slightly with the degree of polydispersity (fig. IT7I) . This 
is surprising since it increases less (maximum 1.16) than 
in the case of the bidisperse media having similar frac- 
tion of fines and size range. We assume that the beads 
are vertically sorted in the flow: the smallest beads at 
the bottom, and progressively larger particles are located 



FIG. 16: Flow and deposit morphologies for an increasing 
polydispersity degree, from top to bottom: batches 2, 3, 4, 
and 8 (see table I for size distributions) (flux: 8.15 g/s, rough- 
ness: 425-600 /im, slope: 25°). 
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FIG. 17: (top) (O) The width and (♦) the velocity of the flow 
seem to be independent of the polydispersity degree. The 
variability of the flows is stronger than errors coming from 
the measurement method (error bars). The vertical axes have 
been chosen to compare these date to the variations induced 
by the range of flux we previously studied (fig. 1131 . (bottom) 
Heights of (0) levees and (•) flow increase with polydispersity 
degree, when (a) the thickness of the channel decreases. X- 
axis is arbitrary: polydispersity increases from left to right 
(X values correspond to batch numbers in table I) (flux: 8.15 
g/s, roughness: 425-600 fim, slope: 25°). 
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higher in the layer, the largest ones being at the top. In 
that hypothesis, the basal smoothing process looses effi- 
ciency when not applied directly to the largest particles, 
as for the case of a bidisperse media. This velocity varia- 
tion has to be compared to the velocity variation induced 
by the flux variation. The range of fluxes and degrees of 
polydispersity used in our experiments are valuable to 
be compared because they both induce a similar varia- 
tion of height. The influence of the flux (fig. I13fl and of 
the polydispersity degree (fig. IT7I) allows a comparison 
of the sensitivity of the three variables (h, v, L) with 
the two parameters (flux, polydispersity). Compared to 
flux-induced variation, the influence of polydispersity is 
negligible and velocity could be considered as constant. 
In the same way, the width of the flow is constant (vary- 
ing less than its fluctuations) with the polydispersity, and 
depends only on the flux, i.e. on the dynamics. In conclu- 
sion, width and velocity vary only slightly with polydis- 
persity, and might be used to establish a relation which 
is independent of the polydispersity. 

The width L seems to be more promising than hi ev ^ e 
for making a link between the morphology of the deposit 
and the dynamics of the flow, which can be transposed 
further to the natural case. However, we will have to look 
for the horizontal length scale of the system, which will 
be the relevant parameter to normalize L, in the same 
manner that h can be normalized by h s t op - 

IV. CONCLUSIONS 

Good accordance between field observations and our 
experiments shows that some small volume pyroclastic 
flows behave as dense granular flows. For these dense 
granular flows down rough inclines, several flow regimes 
have been shown to occur depending on the inclination of 
the plane and on the flux imposed by the source. In each 
regime, the flow obeys different dynamics and is linked to 
a particular deposit morphology. The qualitative obser- 
vation of the deposit allows the determination of the flow 
regime. The determination of the regime fields for a nat- 



ural system could even be a step towards the establish- 
ment of an equivalence between volcanic particles flow- 
ing on a natural substratum and some laboratory gran- 
ular flow conditions. In fact, one limits between regimes 
corresponds to the h s top(9) curve used in laboratory to 
normalize the equations of dense granular flow dynamics. 

In the steady state frictional regime, the flow propa- 
gates with a constant height, width, and mean velocity. 
The associated deposit displays the levee/channel mor- 
phology often observed in the field for unconfined small 
pyroclastic flows. This morphology is without doubt the 
signature of the steady state frictional regime. This mor- 
phology results from static lateral zones and a partial 
emptying of the central part of the flow. The formation of 
the morphology shows that the height of levees is linked 
to the thickness of the flow, and it is possible to quan- 
titatively relate the flow velocity and the morphological 
parameters of the deposit. The method is very promis- 
ing for future modeling of dense pyroclastic flows because 
quantitative information about the flow can be extracted 
only from deposit measurements. A relation has been 
established in the case of the experimental data, for the 
steady state regime associated to the levee/channel mor- 
phology. The forthcoming step could be considered in 
two different ways. First, the relation established in the 
laboratory for a small range of particle sizes could be 
extended to the case of a highly polydisperse material. 
Second, width and velocity could be directly linked, be- 
cause they seem to depend on dynamics of the flow, and 
not on polydispersity. However, we need a relevant pa- 
rameter to normalize the width in order to extend the 
relation to the case of pyroclastic flows. The experimen- 
tal study of the dynamics of granular flows and of the 
formation of the associated deposits morphology could 
be a way to indirectly infer data about the dynamics of 
small pyroclastic flows. 
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